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Abstract 

In this work we present the results of imaging simulations performed with 
the help of the GEANT4 package for the protoMIRAX hard X-ray balloon 
experiment. The instrumental background was simulated taking into ac¬ 
count the various radiation components and their angular dependence, as 
well as a detailed mass model of the experiment. We modeled the meridian 
transits of the Crab Nebula and the Galatic Centre region during balloon 
flights in Brazil (~ —23° of latitude and an altitude of ~ 40 km) and intro¬ 
duced the correspondent spectra as inputs to the imaging simulations. We 
present images of the Crab and of three sources in the Galactic Centre region: 
IE 1740.7—2942, GRS 1758—258 and GX1-1-4. The results show that the pro¬ 
toMIRAX experiment is capable of making spectral and timing observations 
of bright hard X-ray sources as well as important imaging demonstrations 
that will contribute to the design of the MIRAX satellite mission. 
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1. Introduction 


X and low-energy (< 1 MeV) 7 -ray photons coming from astrophysical 
sonrces interact in the Earth’s atmosphere throngh photoelectric effect or 
Compton scattering, never reaching the snrface. However, at nsnal strato¬ 
spheric balloon altitndes (~40km), a major fraction of the X-ray photons 
with energy greater than ~30keV is able to pass throngh the remaining 
external atmosphere. This allows the detection and observation of cosmic 
sonrces above these energies with instrnments carried by balloons. In ad¬ 
dition to photons coming directly from the astrophysical sonrces of inter¬ 
est, radiation from the atmosphere and from other cosmic sonrces (either 
point-like or diffnse) in the held of view of an X-ray telescope interacts with 
the instrnment and the balloon platform, generating an intense backgronnd. 
Therefore, detailed predictions of the hnxes and spectra of cosmic sonrces 
are crncially dependent npon an accnrate knowledge of this backgronnd. At 
balloon altitndes, this is a non-trivial task that has to take into acconnt 
the several particle and photon helds incident npon the experiment and all 
the interactions that occnr in the detectors and snrronnding materials, as 
well as a detailed mass model of the instrnment. By knowing the particle 
kinds, both primary and secondary, their energy spectra, and by following 
the secondary emission they prodnce when they interact, we can simnlate 
the energy depositions on the detectors and identify the contribntion of each 
backgronnd component. This is essential to estimate the instrnment sensi¬ 
tivity as a fnnction of energy. 

In this work, we hrst show the resnits of backgronnd caicniations nec¬ 
essary to make realistic simnlations of the astrophysical observations to be 
carried ont by the protoMIRAX experiment. protoMIRAX consists of a wide 
held coded-mask hard X-ray imager that will operate in a stabilized balloon 
gondola with hne pointing capability. Besides being a prototype designed to 
test several snbsystems of the MIRAX satellite experiment (Braga, J. et al.| 


2004 Braga, 2006) in a near-space environment, protoMIRAX is also capable 


to carry ont spectral observations of bright X-ray sonrces. 

We then present imaging simnlations of known hard X-ray point sonrces 
in selected sky helds to be observed by the experiment in the 30-200 keV 
energy range. We nse standard cross-correlation techniqnes to reconstrnct the 
images with coded masks. Even thongh coded apertnre imaging is inherently 
a low signal-to-noise-ratio (SNR) techniqne, since the whole detector area 
is nsed to measnre both sonrce and backgronnd, it remains the most nsed 
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imaging technique at hard X-ray and low-energy 7 -ray energies due to its 
relative simplicity of implementation and the capacity to image large helds 
of view with reasonable angular resolution. Although the NuSTAR mission 
has extended X-ray optics up to ~ 80 keV (Harrison, F. A. et ah, 2013 


Hailey, C. J. et ah, 2010), coded masks are still important for hard X-ray 


and low-energy 7 -ray wide held imaging instruments. 

Even though protoMIRAX is chiehy a prototype designed to test several 
subsystems of the MIRAX satellite experiment in a near-space environment, 
we show in this work that the experiment is also capable of making spectral 
observations of selected bright hard X-ray sources and perform important 
imaging demonstrations that will contribute to the design of the MIRAX 
mission. 

In section]^ we describe the protoMIRAX experiment, which will be as¬ 
sembled in a gondola to hy on stratospheric balloons. In section we make 
a brief description of the nature of the hard X-ray point sources that we 
have used in our simulations. In section]^ we describe the procedure we took 
to simulate the instrumental background using GEANT4. In section we 
show the calculations we made to simulate the detected fluxes of the sources 
during their meridian passages. In section we show the simulated images 
and discuss their properties. Finally, in section we present the conclusions. 


2. The protoMIRAX experiment 

The protoMIRAX balloon gondola houses the X-ray camera, which is the 
imaging unity of the experiment, and other subsystems. Table [T] presents an 
overview of the experiment and its baseline parameters. 

The camera includes a coded mask based on an extended (4x4) pattern 


of 13 X 13 Modihed Uniformly Redundant Arrays (MURA) (Gottesman and 


Fenimore, 1989), which is placed 650 mm away from a position-sensitive de¬ 


tector plane. The mask elements are made of lead with 20 mm x 20 mm in 
area and 1 mm thickness. Such a design will allow a ~ 1°45' geometrical 
angular resolution and a ~ 20° x 20° fully-coded held of view (FOV). 

The detector plane comprises 169 square (10 mm x 10 mm) Gadmium Zinc 
Telluride (GZT) detectors, 2-mm thick, providing a ~ 30-200 keV optimum 
energy coverage. The lower threshold is dehned by atmospheric absorption 
at balloon altitudes, whereas the upper limit is due to the detector thickness. 
We have developed an electronics acquisition system for the GZTs and our 
current hgure for energy resolution is ~8%@60keV. 
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Table 1: Overview of the ProtoMIRAX balloon mission. 


Mission 

Total mass 

~ 600 kg 

Altitude 

~ 42 km at —23° latitude 

Geometrical area 

169 cm^ 

Effective area @ 50 keV 

50 cm^ (through mask) 

Detector system 

Single detector dimensions 

10 mmxl0mmx2mm 

Material 

Gadmium Zinc Telluride (GZT) 

Number of detectors 

169 (13 X 13) 

Gap between detectors 

10 mm 

Energy range 

30 - 200 keV 

Time resolution 

10 /is 

Coded mask 

Material 

Lead 

Basic pattern 

13 X 13 MURA 

Extended pattern 

2 x2 basic (minus 1 line and 1 column) 

Element size 

20mm X 20mm x 1mm 

Open fraction 

0.497 

Total mask dimensions 

500 X 500 mm^ 

Position 

650 mm from detector plane 

Imaging parameters 

Angular resolution 

1°45' 

Total (fully-coded) FOV 

21° X 21° (FWHM = 14.1° x 14.1°) 


The capabilities of protoMIRAX in the time domain depend upon the 
strength of the source signals, which in turn depends on the detector area. 
The electronic time resolution is 10 ps. For strong sources, we will be able 
to detect time variability at a level of tens of seconds. According to the 
GEANT4 simulations, the expected count rate for the detector plane is ap¬ 
proximately 50 counts/s. At those rates we will have a very low dead time 
and no pile-up problems. 

The balloon gondola has an attitude control system that provides point¬ 
ing capability both in azimuth and elevation with an accuracy of a few arc 
minutes. Fig. shows a computer design of the protoMIRAX experiment. 
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Figure 1: Design of the protoMIRAX payload, consisting of: (1) the pivot with a de¬ 
coupling mechanism for momentum dump; (2) the reaction wheel; (3) the elevation step- 
motor-controlled actuator; (4) star sensors; (5) the MURA-based coded-mask; and (6) the 
electronics bay. The gondola is ~2.4m high, 1.4 m x 1.4 m in cross section, and its total 
mass is around 600 kg 


3. The sample of simulated point sources 

3.1. Crab Nebula 

The well-known Crab Nebula (Ml; a = 5h 34.5m; 6 = -|-22°01’) is the 
remnant of a bright Galactic supernova recorded by Chinese astronomers in 
1054. It is at a distance of 2.2 kpc from Earth and today has an apparent 
magnitude = 8.6. The system encompasses a plethora of regions with 
different structures and emission mechanisms. At the centre of the nebula 
lies the Crab Pulsar, the only neutron star that emits coherent pulses in 
phase all the way from radio to y-rays, with a period of P = 33 ms and P = 
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4.21 X 10 The correspondent spin-down luminosity is Lgd = ~ 

5 X 10^® erg s“^, where I is the moment of inertia of the neutron star (Frank 


et ah, 2002). Only a small fraction of Lgd goes into pulsed emission; the 


majority is carried away by some combination of magnetic dipole radiation 
and an ultrarelativistic wind. This is the source of the energy that powers the 
Crab. It is currently believed that a Pulsar Wind Nebula (PWN) efficiently 
converts the energy of the shocked pulsar wind into the observed synchrotron 
emission, which accounts for ~26% of the total Lsd- For a comprehensive 


review on the Crab, see Hester (2008). 


In X-ray astronomy, the Crab has been considered as a standard candle 
due to its large and nearly constant luminosity. In the literature, the fluxes of 
X-ray sources in a given energy band are often measured in “Crabs”, meaning 
the flux of the Crab in the same energy band. In the range of 2 to lOkeV, 


1 Crab = 2.4 x 10 ® erg cm ^ s ^ 


At X-ray energies above ~30keV, the Crab is generally the strongest 
persistent source in the sky and has a diameter of ~1 arcmin. 


The Crab has experienced some recent flares in the GeV range (Tavani, 


M. et ah, 2011) which were simultaneous or just previous to a particularly 
active period in hard X-rays (Wilson-Hodge, C. A. et ah, 2011). However, 


those surprising variations are of the order a few percent on the 50-100 keV 
energy range and do not make any signihcant difference in the estimations 
made in this work. We will observe the Crab Nebula for flux calibrations 
and imaging demonstrations. 

3.2. Galactic Centre sources 
3.2.1. IEI 740 . 7-2942 

IE 1740.7—2942 is the brightest and hardest persistent X-ray source within 
a few degrees of the Galactic Gentre (GG; a = 17h 45m 40.04s; 6 = 

—29°00'28.1”) and, due to a similar hard X-ray spectrum and comparable 


luminosity to Gygnus X-1 ( 

Liang and Nolan 

1984 

), is classihed as a black 

hole candidate ( 

Sunyaev et al. 

1991 

). Because of the two-sided radio jets as- 


sociated to the source, IE 1740.7- 


■2942 was dubbed the hrst “microquasar” 


(Mirabel et al., 1992; Del Santo et ah, 2007), an X-ray binary whose be¬ 


havior mimics quasars on a much smaller scale. Since the GG direction has 
extremely high extinction, with a hydrogen column density around lO^^cm”^ 
(Sheth et ah, 1996| , a counterpart has not yet been identihed despite deep 
searches made in both the optical and IR bands (Prince et ah, 1991). Ad¬ 


ditional absorption can be attributed to the material of a molecular cloud 
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in which the source is likely embedded (Bally and Leventhal, 1991). The 


source w as studied recently by our group from soft to hard X-rays up to 
200keV (Castro Avila, 2011), showing spectra that can be very well mod¬ 
eled by thermal Comptonization of soft X-ray photons. We will image the 
Galactic Centre region with protoMIRAX with IE 1740.7—2942 in the centre 
of the field-of-view for imaging demonstration and spectral measurements. 

3.2.2. GRS1758-258 

GRS 1758—258, also a microquasar, is one of the two brightest X-ray 


sources near the GC at energies greater than 50 keV (Sunyaev, R. et ah 


1991). Its hard X-ray spectra and variability are similar to that of CygX- 
1 (Kuznetsov, S. et ah, 1996). GRS 1758—258 also does not have yet an 


identified counterpart, so its mass and orbital periods are still unknown. 
The behavior of GRS 1758—258 during the transition between the hard and 
soft spectral states is markedly different from that of Cyg X-1: there is 
an observed time delay of ~ 1 month between changes in its luminosity 
and spectral hardness. This suggests that GRS 1758—258 has two separate 
accretion flows, a thin disk and a halo, which indicates the presence of a 


low-mass companion powered by Roche lobe overflow (Main et ah, 1999). 
GRS 1758—258 will be in the field-of-view of our GC observations, so we have 
included its spectrum here for our imaging simulations with protoMIRAX. 

3.2.3. GXl+4 

The neutron star low-mass X-ray binary system GX1-1-4, the best studied 
accreting pulsar around the GC, is the prototype of the small but growing 
subclass of accreting X-ray pulsars called Symbiotic X-ray Binaries (SyXB), 
by analogy with symbiotic stars, in which a white dwarf accretes from the 
wind of a M-type giant companion. It is a persistent source, but with strong, 
irregular flux variations on various timescales and extended high/low states. 
The system was discovered in 1970 by a balloon observation at energies 


above 15 keV showing pulsations with a period of about two minutes (Lewin 


et ah, 1971). Another balloon observation by Jayanthi et ah (1987) discovered 


pulsations at hundreds of keV. Its peculiar spin history has been the object of 


intensive study (Nagase, 1989), and Makishima, K. et ah (1988) have shown a 


clear transition from spin-up to spin-down behavior at the equilibrium period, 
characterizing a torque reversal episode. Optical flickering with timescales 


of minutes were discovered by Braga et ah (1993) in V2116 0ph, the optical 


counterpart of GX 1-1-4, and optical pulsations at the X-ray frequency were 
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discovered by Jablonski et al. (1997). Pereira et al. (1999) have discovered 


an X-ray frequency modulation of ~304 days in the BATSE/CGRO data, 


whereas infrared spectroscopy by Hinkle et al. (2006) have shown a single- 
line spectroscopic binary orbit with a period of 1161 days,. If this is the 
orbital period, this system is the widest known low-mass X-ray binary. More 
recent studies (Gonzalez-Galan et al. [2012 ) have used data from BeppoSAK, 
INTEGRAL, Fermi and Swift/BPCT to show that the source continues its 
spin-down trend with a constant change in frequency, and the pulse period 
has increased by ~50% over the last three decades. 

We will have GX1-1-4 in the field-of-view of the GG protoMIRAX obser¬ 
vations. The detection of the source will help to demonstrate the imaging 
capabilities of the telescope and make an updated measurement of the spin 
period of this very interesting and peculiar object. 


4. Simulation of the background at balloon altitudes 

To compute the protoMIRAX instrumental background at balloon alti¬ 
tudes we made use of the well-known radiation-interaction code GEANT4, 


developed by GERN (Allison, J. et ah, 2006 Agostinelli, S. et ah, 2003). We 


have built a procedure to take into account the main particles (photons, pro¬ 
tons, electrons and neutrons) coming from all directions in the atmospheric 
environment at balloon altitudes (~ 40 km). 

In the case of the external electromagnetic radiation, we considered the 
anisotropic atmospheric X and y-ray spectra and the cosmic diffuse contri¬ 
bution coming through the telescope aperture. We have divided the contri¬ 
butions in 4 zenith angle {z) ranges from 0° to 180°, in the 0.024 — 10 MeV 
energy range. Table shows the different power-law spectra adopted in each 
angular range. 

The proton spectrum consists of two components (see Table [^. The 
Erst one (primary) is the most energetic and corresponds to the cosmic rays. 
They will come only from the upper hemisphere, since they are from cosmic 
origin. The second component (secondary), less energetic, is the product of 
the scattering of the primary protons in the atmosphere. Therefore, these 
protons can come from any direction. We calculated the primary proton 
distribution in the range 10^ —1.5 x 10^ MeV, while for the secondary protons 
we considered the range 5 — 6600 MeV. 

The electrons appear as a consequence of the interaction of cosmic rays 
in the atmosphere. Their spectrum in the 1 — 10 MeV energy range is given 




















Table 2: Incident photon spectra in units of photons cm ^ s ^ MeV ^ in the 0.024 — 10 
MeV energy range, for zenith angles between 0° and 180°. Taken from Gehrels (1985). 


Spectrum 

Angular range 

r] 

Energy range 
[MeV] 

2.19 X lO'^E^-'^ 

0 - 65 

0.024 - 0.035 

5.16 X 

0 -65 

0.035 - 10 

0.085 

65 -95 

1 - 10 

0.140 

95 - 130 

1 - 10 

OMOE-^-^^ 

130 - 180 

1 - 10 


Table 3: Primary and secondary proton spectra as a function of energy. Taken from 


Dean et al. 

(2003 

1 

Spectrum 

[protons cm“^ s“^ MeV“^] 

Energy range 
[MeV] 


Primary protons 


1.3 X lO^E-2 

10^ - 3 X 10^ 

1.8 X 10®E-2-^ 

3 X 10^ - 1.5 X 10^ 


Secondary protons 


1.475 X 10-®E2 

5 - 16 

3.78 X IQ-^E^ 

16- 100 

6.43 X 10-3^-0-61 

100 - 300 

1.78 E-i-6 

300 - 6600 


Table 4: Neutron spectra as a function of energy. Taken from Gehrels (19851 

Spectrum 

[neutrons cm-^ s-i MeV-i] 

Energy range 
[MeV] 

0.26E-1-3 

6 X 10-2^-0-65 

0.8 - 10 

10 - 100 


by: 


dN 

dE 


= 1.4 X IQ-^E 


-2 77 - 1.8 


cm ^ sr ^ s ^ MeV ^ 


Table shows the neutron spectrum as a function of energy (Gehrels 


1985). Neutrons are secondary particles, since they appear as a consequence 


of the interaction of high energy protons and alpha particles in the atmo¬ 


sphere (Armstrong et ah, 1973). 
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In order to carry out the background simulations using GEANT4, we have 
built a detailed mass model of the protoMIRAX experiment including the 
detector plane, the collimator, the coded mask, the support structure and 
the shielding materials. This model, together with all photon and particle 
helds described above, were used as inputs to the program. 

GEANT4 takes into account all the physical processes by which particles 
interact with matter. The user can choose which processes are going to be 
considered in a given simulation. Those processes include pair production, 
photoelectric effect, Gompton scattering and Rayleigh scattering for photons, 
multiple scattering, ionization, Bremsstrahlung and Goulomb scattering for 
electrons, and elastic and inelastic processes and Bertini cascades for protons. 
The processes considered for neutrons include those for protons, plus neutron 
capture and hssion. 

Fig. I^shows the spectrum of each background particle species considered: 
photons, electrons, protons and neutrons. By adding all the interactions that 
deposit energy onto the protoMIRAX detector plane by direct or indirect 
hits, considering all the photon and particle helds used in this simulation, we 
can produce the total background spectrum expected for the protoMIRAX 
balloon hights. This is shown in Fig. for an integration of 8h at an 
atmospheric depth of 3 g cm“^ at a latitude of —23° over Brazil. 

The total number of simulated background counts is 1,453,260, which 
corresponds to 50.46 ± 0.04 counts/s. The background drops considerably 
above 40 keV and stays reasonably hat at a level of ~ 1.6 x 10“^ counts"^ 
s“^ keV“^ up to 200 keV. In Fig. we show the spatial distribution of 
the background across the detector plane. The background is reasonably 
smooth, with slightly higher counting rates near the edges and corners due 
to geometrical factors. The details of the background simulations, including 
the relative importance of each input spectrum and the specihc properties of 
the GEANT4 input hies we have used {e.g. the mass model), are the subject 
of another work in preparation that will appear elsewhere soon. 

5. Simulations of the incident fluxes from the point sources 

The X-ray/low-energy y-ray photon spectrum of the Grab Nebula is 
known to be well described by a power-law 

Fq = AqE~'^ photons cm“^s“^keV“^, (1) 
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Figure 2: Background spectra estimated for each particle species (electrons, photons, pro¬ 
tons and neutrons) for protoMIRAX based on GEANT4 simulations at balloon altitudes 
at a latitude of —23° S. The error bars are la statistical uncertainties. We have taken into 
account the activation of the detector materials in all cases. Six lead fluorescence lines are 
visible in the photon and electron spectra at 72.8, 74.9, 84.4, 84.9, 87.3 and 88 keV (values 
taken from http://www.kayelaby.npl.co.uk/atomic_and_nuclear_physics/4_2/4_2_l.html) 


with Aq = 14.4 and 7 = 2.2 (Sizun, P. et ah, 2004). The fluxes at each 


energy are attenuated due to absorption in the atmosphere, so they vary as 


F{E) = Fo{E) 


( 2 ) 


where fi/p{E) is the absorption coefficient of the atmospher^at energy E, 
given in cm^ g“^, x = 2.7 gcm“^ is the atmospheric depth at the expected 
balloon altitude (~ 40 km), and 2 ; is the zenith angle. 

In order to calculate the total number of counts from the Crab in the 
protoMIRAX energy range during its meridian passage, we have to integrate 
over its path in the sky. We considered a 4-hour period divided equally before 
and after the meridian transit, since the source would be too low in the sky 
before and after that. 


^http://physics.nist.gov/PhysRefData/XrayMassCoef/ComTab/air.html 
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Figure 3: Background spectral contributions estimated for protoMIRAX based on 
GEANT4 simulations at balloon altitudes at a latitude of —23° S. Red corresponds to 
photons, blue to protons, cyan to neutrons, pink to electrons and black represents the sum 
of all the contributions. The error bars are la statistical uncertainties. 


First, we consider the dependence of the zenith angle with time. Let 0 
be the geographic latitnde (</> = —22°39' at Cachoeira Panlista, Sao Panlo, 
Brazil), 6 = +22°01' the declination of the sonrce, a = 05h34.5m the right 
ascension and H = t — a the honr-angle, t being the sidereal time (LST). 
The zenith angle 2 ; is given by 


cos z{t) = sin 0sin 5 + cos (j) cos 6 cos {t — a). 


( 3 ) 


Given the expression for z{t), we can compnte the nnmber of photons N 
coming from the Crab reaching the detector of area S' in a time interval T 
in the detector energy range Em, 


Emax- It will be 


N = 


AqE '*'exp 


■ —(E) xsecz 
P 


dEdtdS. 


( 4 ) 


For the protoMIRAX experiment, the total energy range is 30 to 200 keV 
and the area to be considered for this calcnlation is ~ 84 cm^, since this 
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Figure 4: Count distribution of the total background over each pixel on the detection 
plane in the detector energy range 30 — 200 keV. The observation time was 8 h at balloon 
altitudes. Each pixel corresponds to one planar CZT detector (1 cm^). 

corresponds to the number of detectors that are not covered by the closed 
elements of the mask during an observation of a point source. This is the 
case because, for the 13 x 13-element MURA pattern used (see section 5), 85 
cells are closed and 84 are open. For the Crab simulation, we considered 
it to be at the centre of the £eld-of-view, so there are no losses due to the 
collimator response. 

We performed the integration numerically by using sufficiently small en¬ 
ergy and time bins so as to get a good approximation, given the variations of 
the spectrum with energy and of 2 ; with time. After calculating the number 
of counts N for each energy bin, we can run a GEANT4 simulation by shin¬ 
ing the corresponding absorbed Crab spectrum over the protoMIRAX mass 
model. 

Fig. (left) shows a map of the total counts coming from the Crab in 
each detector for the 4-hour meridian passage simulation. This Crab map 
is actually a “shadow” of the mask projected onto the detector plane, since 
the flux is spatially modulated by the mask (see section 5). Each pixel of 
the map corresponds to one detector. The total number of counts due to the 
Crab is 103331. When we add a 4-hour simulated background, we get the 
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complete shadowgram shown in Figj^ (right). 


Crab Crab + Background 



Figure 5: Crab (left) and total (Crab + background) (right) count distribution over the 
detection plane in the detector energy range 30 — 200 keV. The observation time was 4 h. 

Following the same procedure as the one used for the Crab Nebula, we 
considered the emission from three sources near the GC: IE 1740.7—2942, 
GXl+4 and GRS 1758—258, as described in section]^ The parameters of 
the photon spectra of each of these sources are presented in Table We 
have used these spectra and the zenith angle variation of the GG direction 
to perform the integral of equation (|^ for each source. 

Table 5: Photon spectra of the source sample near the GC. The spectra are in units of 
photons cm“^ s“^ keV“^ 


Source 

Spectrum 


Reference 


IE 1740.7-2942 

iq-4 E -1-^5 

lOOkeV 

Grebenev et ah 

(1995 

) 

GRS 1758-258 

4-6 X 10-5 

Sunyaev.et ah (1991b) 

GXl+4 

5-1 X 10-\,,f,^ 


Dieters et ah ( 

199ll 



In this case, since the GG transits will occur at < 6° from the zenith 
at the latitudes of the protoMIRAX programmed balloon flights, we have 
considered integration times of 4 hours before and after the transits. In or¬ 
der to maximize the signal-to-noise ratio (SNR) of IE 1740.7—2942, we have 
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centralized in the image in this source. Fig. shows the correspondent 
shadowgrams for the GC simulation. It is interesting to note that the shad¬ 
owgram, in this case, bears almost no resemblance to the mask pattern (see 
section 5), since the overlapping shadowgrams from the three point sources 
interfere in a complex way. 



Figure 6: The three sources we studied in the GC (left) and the total counts distribution 
(3 sources + background) (right) over the detection plane in the detector energy range 
30 — 200 keV. The observation time was 8 h. 

It is noteworthy that in the actual observations we will have to take into 
account the sky rotation in the £eld-of-view during the integration times, 
due to Earth’s rotation. For this simulations, we assumed that this effect will 
be properly dealt with by rotating the shadowgrams obtained during short 
integration times, during which the source positions will change by a small 
angle compared to the instrument resolution, and the rotated shadowgrams 
will be stacked and added up with the proper orientations. 


6. Image reconstruction 

In a coded-mask experiment, the mask is composed of elements which 
are transparent (open elements) or opaque (closed elements) to radiation. 
Photons hitting an open element will pass through the mask and reach a 
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position-sensitive detector plane, and photons hitting a closed element will 
be absorbed. In this way, a pattern of mask “shadows” is projected on the de¬ 
tector plane, which actually corresponds to the superposition of many images 
of each object (one image for each open element). In order to determine the 
intensity and position of each observed object in the held of view, a decoding 
and reprocessing of the shadowgram is needed through the implementation 
of numerical algorithms (Caroli et ah, 1987). This can be a very complex 


process in the general case of extended sources at different distances from the 
instrument. In the present case of point sources at inhnity, the reconstruction 
procedure is somewhat simplihed. Essentially, the reconstructed image is the 
cross-correlation between the detector count distribution map and a decoding 
function that mimics the mask pattern (Fenimore and Cannon, 1978). 

The coded mask used by protoMIRAX is shown in Fig. It is a cyclic 
(4 X 4) repetition of a 13 x 13 MURA pattern ( Gottesman and Fenimoi^ 
1989), minus one line and one column. By using a MURA pattern we get 
“perfect” imaging (no intrinsic noise) and the incomplete repetition of the 
basic pattern in all directions ensures that we have no ambiguities in the 
locations of the sources in the field-of-view, since each object always projects 
a distinct shadow on the detector plane. Since the collimators completely 
block X rays coming from directions not covered by the mask, we have no 
partially coded field-of-view in this experiment. 

Based on Poisson statistics, the signal-to-noise ratio (SNR) of a peak 
{i,j) in a reconstructed coded-mask image, which determines the sensitivity 
of these systems, can be shown (Gottesman and Fenimore, 1989) to be given 
by 


SNRij 


Ns 

VNs + Nt ’ 


( 5 ) 


where Ng is the number of net source counts (corresponding to source photons 
that came through the mask openings and hit the detectors), and Nt is the 
total number of detector counts, both for the same integration time. 

We show in Fig. [^the reconstruction of a simulated image of the Grab 
Nebula region produced by cross-correlating the shadowgram of Figj^ (right) 
with the mask pattern. The total £eld-of-view of the image is 21.7° by 21.7° 
and the Grab is detected in one sky-bin of 1°45'. By dividing the peak value 
by the standard deviation of the off-peak region, we get a SNR of 109, which 
is 91% of the theoretical value given by equation 5. This result shows the 
good performance of the imaging system of the protoMIRAX experiment. 
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Figure 7: The coded mask pattern used by protoMIRAX. It is a cyclic (4 x 4) repetition 
of a 13 X 13 MURA pattern, minus one line and one column. The basic cells are made of 
Imm-thick lead and measure 10 mm x 10mm. 

Also, our results indicate that we will be able to make a 5a detection of ~1- 
Crab sources in 30 s, which will be important to study high-resolution time 
variability of transient sources. 

For the Galactic Centre region simulation, we considered IE 1740.7—2942 
to be in the centre of the £eld-of-view. The angular distance from this source 
to GX 1-1-4 is 0 = 5.65°, and to GRS 1758—258 is 0 = 5.53°. These sep¬ 
arations are signihcantly greater than the protoMIRAX angular resolution, 
which shows that the instrument is capable of measuring hard X-ray spectra 
and light curves in relatively crowded fields like the GC. 

Fig. 1^ shows the reconstructed image of the GC region. To calculate the 
SNR of each pixel, we hrst made the identihcation of the peaks by selecting 
pixels above 3 standard deviations of the entire image, after subtracting the 
DC level. We then calculated the SNRs by dividing the peak values by the 
standard deviation of the off-peaks region. The detection levels are approxi¬ 
mately 60% of the purely statistical values. The signal loss is a combination 
of 3 effects: non-uniformity of the background across the detector plane (see 
Fig. 3), the collimator response and the existence of gaps between detectors. 
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Figure 8: Simulated image of the Crab region as seen by protoMIRAX for 4 hours at 
balloon altitudes. The Crab is detected at a 109cr level. 

7. Conclusions 

In this work we presented results of simulated images of regions with 
bright hard X-ray sources with the protoMIRAX balloon experiment. pro¬ 
toMIRAX is a prototype designed to be a scientific and technological pathfinder 
for the satellite wide-held hard X-ray imager MIRAX. 

We introduced the results of GEANT4 simulations of the physical con¬ 
ditions into which the proto-MIRAX experiment will be subject to in its 
stratospheric balloon hights, in order to predict the performance of the tele¬ 
scope during the observation of bright cosmic sources in hard X-rays. We 
modeled the instrumental background expected during hights at altitudes of 
~ 40 km over Southeast Brazil (~ —23° latitude) by simulating the inter¬ 
action of incident radiation and particles (photons, electrons, protons and 
neutrons), cosmic and atmospheric, with the experiment materials. 

We then produced images of two regions which will be observed by pro¬ 
toMIRAX in its hrst balloon hights for calibrations and imaging demonstra¬ 
tions: the Crab Nebula and a Galactic Centre held with 3 bright sources: 
IE 1740.7—2942, GRS 1758—258 and GX1-1-4. The expected huxes of the 
sources at balloon altitudes, in several energy bands, were added to the 
background in order to produce shadowgrams of event distributions over the 
detector plane. 
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Figure 9: Simulated image of the Galactic Centre region as seen by protoMIRAX for 8 
hours. The SNRs for IE 1740.7-2942, GRS 1758-258 and GXl+4 are, respectively, 26, 
10 and 9. We can see that GXl+4 counts are distributed in two sky bins. 

We simulated the Crab Nebula image considering a meridian passage of 4 
hours, during which the highest source elevation will be 45°. We have shown 
that we will be able to reach a very good SNR of 109, which will make the 
Crab a reliable calibration source for the experiment. 

In the case of the GC, the three sources lie near the centre of the field 
of view and were imaged together, with IE 1740.7—2942 in the centre. The 
results show that we can make observations of crowded fields like the GC 
with protoMIRAX and provide flux and spectral information without source 
confusion. 

The protoMIRAX mission, as other scientihc balloon missions, plays a 
very important role in testing new detector technology and imaging systems 
in a near-space environment. As a pathhnder experiment for the MIRAX 
satellite mission, protoMIRAX is an invaluable source of information for 
technology development and demonstration. 
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